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A quantum-kinetic equation for the description of 
internal dynamics of multilevel atomic systems 
moving through a target matter 

Voskresenskaya 1 

Joint Institute for Nuclear Research, Dubna, Moscow Region, 141980 Russia 

Abstract. The quantum mechanical consideration of a passage of relativistic 
'^T)' 1 elementary atoms (EA) through a target matter is given. A quantum- kinetic 

•O equation for the density matrix describing their internal state evolution at EA 

| rest frame is derived. 

1 Introduction 

! For the interpretation of the data of DIRAC experiment [1, 2] which aims to measure 

the lifetime of hydrogenlike EA consisting of 7r + and tt~ mesons (^2^ atoms) one 
needs to have the accurate theory for the description of internal dynamics of the 
atoms moving through a target matter. 

During their passage through the target (pionium atoms) interacts with 
target atoms that causes the excitation, deexcitation or ionization of the A2-K- To 
describe these variations of A 27T internal states the authors of [3] proposed a set 
of kinetic equations for the probabilities to find the pionium atom in the definite 
^ ■ quantum state at some distance from the point of A 2 tt production. 

It is clear that such "classical" description is approximate because does not take 
into account the possible interference (quantum) effects. These last can be included 
r-| ! in consideration only in the framework of density matrix formalism. 

> : 

2 Derivation of a density matrix kinetic equation 

S ; 

The derivation of a kinetic equation for the density matrix of fast atomic systems 
passing through a target matter can be given at target rest frame [4], but more 
simple these equation can be obtained at rest frame of the EA. 

At this frame the target moves with the velocity Vq and the electromagnetic 
field produced by target atoms is described by 4- vector potential = (Q,A), 
A=(v /c)<f>. 

The scalar potential $ interacts with the charges of mesons and the vector po- 
tential A^ with their currents. Because the typical velocities of the particles forming 
EA are of order ac <C c (a is the fine structure constant), we will neglect the term 
proportional to the current in the Hamiltonian (see [5]). 

Then internal dynamics of relativistic EA (later, for definiteness, "of pionium 
atoms") is described by the Schrodinger equation 
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with the Hamiltonian of the form 

H = H + V{r,t), H = T + V (f) (2) 

and 

T = -A/2// = -{d/dr) 2 /2/i . (3) 

Here, Vo(r) are the potential energy of pion-pion interaction and V{f,t) is the po- 
tential energy of interaction between the pionium and the target atom. 

We will suppose that the positions of atoms inside the target are not varied 
during the interaction of target with the pionium atom (the so-called "frozen" target 
approximation). Then 

V(f) = e£> {nit) - r/2) - $ {nit) + f/2)] , (4) 

i 

fi{t) =n{t ) +v {t-t ), (5) 
${R) = 7 W^ 2 + 7 2 (^) 2 , (6) 

R = f l {t) T r/2, 7 = Wl-^. (7) 

Here, $o is the potential of the target atom at it's rest frame, and we have put the 
origin of the coordinate system to the center-of-mass of pionium. 

Thus, the solution of the Schrodinger equation (1) depends on "frozen" positions 
fi{t ) of the target atoms 

i>{r,t) =4>(f,t; {fi{t )}). 
The density matrix of pionium is defined as follows: 

p{r,r';t) = (^(r-t;^^)})^^,^^^)}))^^, (8) 

where (){ff(*o)} m eans averaging over all possible positions of target atoms. 

Let to be the point of time when moving target meet the pionium atom and 
ip{r,t ) is the value of pionium wave function at this time. Then at t > t 

^(r,t;{fi(* )})= J G(f,f \t,t \{fi(to)}^ipi(fo,to)df , (9) 

where G is the Green function of Eq. (1). 

According to [6], it can be expressed in terms of the path integral 

G{r,r ;t,t ;{f i {t )}) = J Dr{t) exp(iS') , (10) 

with 

S = S + S U (11) 
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S = J dt'L (v(t'),f(t')) , S 1 = -J dt'V(f(t'),t') , (12) 



to to 

L (v(t'),r(t')) = »v 2 (t')/2-V (r(t')), (13) 

v(t') = df(t')/dt' . 
It can be shown (see [9]) that 

Si = -E{x(6< + W/2)-x(6i-^i)/2)}^-*i), (14) 

i 

where 

oo 

e 



x(k) = - / ®{^b 2 ± + z*)dz, (15) 



-oo 



6 ± = 6 i± ^, fi = io+ *4M, (16) 



6» = r*(*o)± = ri(t ) - ^° Tl 2 ^ ■ v , (17) 



s(U) = r(U)± = m - V ° ' n } to) ■ vo , (18) 

the Heavyside step function $(t) is for t < and 1 for £ > 0. 

Substituting (10)-(18) into (8) and performing the averaging over the "frozen" 
positions of target atoms with the help of the prescription of [7, 8], one can get the 
following representation for the density matrix: 

p(f, f t) = J G(r, f r , fg; t, t )ipi(r , t )^*(^o, t Q )df di% , (19) 

G(f, f '■ f , f£ i, t ) = | Df{t)Df'{t) exp(iS - W) , (20) 
t 

S = Jdt'{L (v(t'),r(t>))-L (v'(t'),r'(t'))} , (21) 

to 

t 

W = v Ql n J dt'Sl(s(t')J'(t')), (22) 

to 

n(s(t'),s'(t')) = |rf 2 6{l-exp(^(6,s(t , ),s'(t')))} , (23) 



with 



- X (b + s'(t')/2)+(b-s'(t')/2). (24) 

Here, no is the number of atoms in the unite volume of target at it's rest frame, s 
and s' are the transverse parts of the vectors f and r'. 

From Eqs. (19)-(22) it easily derive (see [6]) the following equation for the density 
matrix: 

i dp(r^ 1 t)_ = Ho(f)p(f y. t )-H (r')p(r,r';t) 

-iv ^n ^(s, s')p(r, r'\ t) , (25) 

where the last operator term describes the Coulomb interaction between EA and the 
target atoms with account of all multiphoton exchanges. Using a generalized optical 
potential of the form V opt (s, s') = kQ(s,s'), where k = —ivojiio, we can represent 
this term as V opt p(t). 

The form of Eq. (25) is similar to the form of Eq. (116) in Ref. [10] describing 
the internal dynamics of multilevel atoms in laser fields, where the last term Tp 
describes the contribution of the spontaneous relaxation. 

The equations of motion for the density matrix elements 

Pik = J ip*(r)^k(r')p(r,r')drdr' (26) 

looks like as follows: 

dpl ^ = iA ik p ik (t) - v -fn ^ik,imPim{t) , (27) 

l,m 

where 



V ikM = J r t (^M^Mr'W m (r'Ms, s')drdf , (28) 
the EA wave functions ip^k) and binding energies e^k) obey the Schrodinger equation 

H ^i(k) = £i(k)^i(k) ■ (29) 
Taking into account the lifetime of the EA, we can obtain 



dpikjt) 
dt 



1 

i(e k - Si) - -(Ti + T k ) 



Pik (t) - Vo'YTlo ^ik,lmPlm if) , (30) 
l,m 



where = l/r^ is the EA levels width (for details see [4]). 

Properties of the solution of (25) will be discussed in the following paper. 
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